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ABSTRACT 

We investigate the relation between star formation (SF) and black hole accretion luminosities, using a sample 
of492 type-2 active galactic nuclei (AGNs) at z < 0.22, which are detected in the far-infrared (FIR) surveys with 
AKARI and Herschel. We adopt FIR luminosities at 90 and 100 pm as SF luminosities, assuming the proposed 
linear proportionality of star formation rate with FIR luminosities. By estimating AGN luminosities from 
[O III] A5007 and [OI] A6300 emission lines, we find a positive linear trend between FIR and AGN luminosities 
over a wide dynamical range. This result appears to be inconsistent with the recent reports that low-luminosity 
AGNs show essentially no correlation between FIR and X-ray luminosities, while the discrepancy is likely due 
to the Malmquist and sample selection biases. By analyzing the spectral energy distribution, we find that pure- 
AGN candidates, of which FIR radiation is thought to be AGN-dominated, show significantly low-SF activities. 
These AGNs hosted by low-SF galaxies are rare in our sample (~ 1%). However, the low fraction of low-SF 
AGN is possibly due to observational limitations since the recent FIR surveys are insufficient to examine the 
population of high-luminosity AGNs hosted by low-SF galaxies. 

Subject headings: galaxies: active — galaxies: star formation — infrared: galaxies 


1. INTRODUCTION 

In decades the connection between galaxy evolution and the 
growth of supermassive black holes (SMBHs) has been one of 
the main topics in extragalactic research. The tight correlation 
of black hole mass, Mbh. with galaxy properties, e.g., stellar 
velocity dispersion, <r* (e.g., Magoman et al.|1998| Marconi| 


& Hun t|2003] [Woo et al. |20 lOf [2013|7 suggests the coevolu 

tion of galaxies and SMBHs although the physical link be¬ 
tween them is yet to be clearly revealed. Various observa¬ 
tional studies have been devoted to investigating the nature 
of the coevolution. For example, the redshift evolution of 
the Mbh-ct* relation, representing a cumulative growth his¬ 
tory, has been investigated mainly using type-1 active galac¬ 
tic nuclei, AGNs (e.g., Woo et al. |2006 2008 Merloni et al. 


|201 0[ [Schramm & Silverm an 2013). The connection between 

on-going star formation (SF) and AGN activity is also one 
of the observational signatures, revealing the connection of 
the growth of stellar mass and the BH growth at the observed 
epoch (e. g., [S anders et al. |1988 Cid Fernandes et al. |2001[ 


Kauffmann et ak[ 2003 |Netzer||2009[ [Alexander & Hickox 

2012| Rosario et a .|2012 1 . 


Various theoretical frameworks have been suggested to ex¬ 
plain the AGN-SF link and also reproduce the Mbh-<z* rela¬ 
tion. For example, based on the smoothed particle hydrody¬ 
namic /V-body simulations of gaseous galaxy mergers, [Blecha 


[et al.| ( |20TT] i showed simultaneous bursts of SF and BH accre¬ 
tion (see also |Hopkins & Quataert||2010| ). Such theoretical 
models predict a positive correlation between SF and AGN 
activity, as consistent with the results of several observational 
studies (e .g., Netzer etal.|2007 jNetzer 2009[|Diamond-Stanic| 


|& Riekel|2012[ |Rovilos et al. |2012| |Karouzos et al.||2014j i7 

Note that other studies reported that there is a time lag be¬ 
tween SF and AGN phases (e.g., IDavies et al.|2007l|Wild et| 
|al.|2010[|Matsuoka et al.|20TT]|Cen|2012[|Hopltins|2012|i. 
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While the AGN-SF link can be investigated in various as¬ 
pects, the direct comparison between AGN and SF luminosi¬ 
ties, i.e., the Lagn-Lsf relation, is the most simple approach. 
Since the SF luminosity corresponds to on-going growth of 
galaxies and the AGN luminosity reveals the current growth 
of SMBHs, the AGN-SF connection can be directly traced at 
the observed epoch. A correlation between SF and AGN lumi¬ 
nosities has been reported in the previous studies, indicating 
that luminous AGNs are hosted by highly star-forming galax¬ 
ies. Based on the combined sample of local t ype-2 AGNs 
and quasars at 0.1 < z < 3, for example, Netzer](20091 found 
there is a good correlation between SF and AGN luminosities, 
albeit with substantial scatters (see also [Netzer et al.||2007[ 
[Lutz et al.|2008[|Woo et al.|2012| i. Recently, Tommasin et al.| 
( |2012| ) have found a correlation between Lagn and Lir of low- 
ionization nuclear emission-line regions (LINERs). These re¬ 
sults suggest that BH activity is connected with SF. 

Based on the deep Herschel imaging of the X-ray sources at 
0.2 < z < 2.5 in the fields of GOODS and COSMOS, |Rosariol 
et ak| (2012| have reported a correlation between AGN and 


FIR luminosities. In their study, luminous X-ray AGNs at 
Z < 1 show a correlation between AGN and FIR (i.e., 60 /jm) 
luminosities, as similarly presented by earlier works (e.g., 
jNetzer et al.||2007] [Lutz et ak|2008[ |Netzer||2009| ) while the 
correlation flattens or disappears at z > 1 (see also, e.g. ,|Lutz| 
et al. 2010 Shao et al. 2010| Hatziminaoglou et al. 2010| 


Harrison et al.||2012| Page et al.| 2012| >. In contrast, they 
claimed that low-luminosity AGNs show essentially no cor¬ 
relation between FIR and AGN luminosities at all redshift. 
The enhanced SF for given AGN luminosity of their low-z 
X-ray AGNs (0.2 < z < 0.5) seems to contrast to the find¬ 
ing of [Netzer] ( |2009j ) that local type-2 AGNs (z < 0.2) show a 
positive correlation between SF and AGN luminosities. This 
discrepancy may be caused by observational biases, e.g., the 
Malmquist and sample selection biases, and from the mea¬ 
surement uncertainties in SF and AGN luminosities. It is also 
possible that for a given AGN luminosity, galaxies with lower 
SF luminosity may be undetected due to the FIR flux limit. 
Therefore, in order to fully reveal the connection between 
AGN and SF activities it is crucial to examine potential bi- 
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ases, which may affect the Lfir-Lagn relation. 

FIR luminosity is often used as a SF indicator since the 
rest-frame FIR emission is mainly from the host galaxy while 
the AGN contribution to FIR, ~ 50-150 pm, is due to the 
Rayleigh-Jeans tail of an AGN-heated dust component (e.g., 
Netzer|2009||Mullaney et alT 2012| Rosario et al. 2012 1 . In the 
case of AGN luminosity. X-ray luminosity can be used with a 
proper bolometric correction. However, X-ray data with suf¬ 
ficient depth is often not available. Instead, emission lines 
from the narrow-line regions (NLRs), e.g., H/3, [O III]A5007, 
[O I]A6300, and [O IV]A25.89//m lines, are often used as 
a proxy for AGN_bolometric luminosity (e.g., |Netzer|l2009 ; 
[Diamond-Stanic & Rieke|2012j l. 

In this paper, we investigate the AGN-SF connection for a 
sample of type-2 AGNs at 0.01 < z < 0.22 selected from the 
Sloan Digital Sky Survey (SDSS), using AKARI and Her- 
schel data. In Section [2] we describe the sample selection 
and the data. Section pfpresents the main results and Sec¬ 
tion]?] provides discussion and interpretation. The summary 
and conclusion are given in Section |5] We adopted a concor¬ 
dance cosmology with (DmMa) = (0.3,0.7) and Hq = 70 km 


Mpc 


2. SAMPLE AND DATA 


In this study, we mainly focus on the local type-2 AGNs 
selected from SDSS, for which FIR data are available. Us¬ 
ing SDSS spectroscopic data and FIR survey data, we inves¬ 
tigate the relation between AGN and FIR luminosities in the 
wide dynamic range. We also collect multi-wavelength data, 
e.g., ultraviolet (UV) and mid-infrared (MIR), to examine the 
characteristics of the galaxies in the sample. In this section, 
we describe our sample selection and multi-wavelength data. 

2.1. Type-2 AGN Sample 

We selected type-2 AGNs at 0.01 < z < 0.22 from the MPA- 
JHU SDSS DR7 galaxy catalog 3 , including 927,552 galax 
ies, based on the BPT diagnostic diagram (e.g., Baldwin et 


al.|1981) using the [O III] A5007/H/3 and [N II]A6584/Ha flux 
ratios with the classification scheme in |Kewley et al.| ( |2006| l. 
Note that we excluded composite objects since they are unre¬ 
liable in estimating AGN luminosities from narrow-emission 
lines due to the contribution from SF. In the selection process, 
we also adopted following criteria: reliable redshift measure¬ 
ment (i.e., ^warning = 0) and high signal-to-noise ratios (S/N) 
for emission lines which were used in AGN classification, i.e., 
S/N > 6.0 for [O III] A5007, Ha, and [N II]A6584 lines, S/N 
>3.0 for the H/3 line. As a result, we obtained 35,945 type-2 
AGNs. 

Note that our sample contains both Seyfert 2s and LINER 
2s, although the physical connection between them is not clear 
as discussed in various studies. Especially understanding the 
ionization process of LINERs is essential in our study because 
it directly relates to the estimate of AGN luminosities (Sec¬ 
tion |3.2) . Mainly there are three considerations for LINER 
2s: (1) hot stars, i.e., post-AGB star and blue stars showing 
similar line-flux ratios to LINERs, (2) shock ionizations in¬ 
stead of photoionozations, and (3) low-ionization parameters 
of LINER 2s. In order to consider these issues, we divide 
our sample into two subsamples, i.e., Seyfert 2s and LINER 
2s (Section |2.2) . To remove stars showing similar line-flux 
ratios to LINERs, we identify the so-called [O I]A6300-weak 

’ http://www.mpa-garching.mpg.de/SDSS/ 
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Fig. 1.— FIR luminosity of the AKARI-detected (90 /im; small circles) 
and Herschel -detected sources (100 //m; stars) as a function of redshift. The 
sample is color-coded based on the luminosity of [O III] A5007 line, as la¬ 
beled at the top left. The dotted and dashed lines represent the 5cr-detection 
limits of the AKARI/FIS and the Herschel !PACS surveys, respectively. 

LINERs, which are considered non-AGNs. By adop t a c ri- 
terion, i.e., [O I]A6300/Ha <1/6 (e.g., Filippenko & Ter-| 
|levich|1992]|Alonso-Herrero et al,|2000) , we exclude LINER 
2s (see Section [L2| . Regarding shock ionizations, we expect 
that shocked gas should lead to higher excitation UV spectra 
than photoionized gas. Thus, strong UV emission lines, e.g., 
C IVA1549, are expected. However, some observational stud¬ 
ies reported featureless UV spectra of LINERs. For example, 
Barth et al ( |1996 1997) presented UV spectra of local LIN¬ 


ERs that high-excitation lines are not detected, meaning the 
fast shock models are poor matched to the observed spectra 
(see also |Maoz et ~aI7||1998 [Nicholson et al,|p 998| |Gabel et| 


|ak||2000[ 


Maoz et al. 

1998 

Sabra et al. 

2003 


In this study, therefore we as- 


jects. Moreover, the low-ionization parameter of LINERs is 
another concern to consider. Since the ionization parameters 
of LINERs are systematically smaller than that of Seyferts, 
there is large uncertainties in determining AGN luminosity 
adopting the same method used for Seyferts. However, if we 
use the method based on both [O III]A5007 and [O I]A6300 
lines, we can correct for the ionization effect of low-ionization 
sources such as LINERs (see Section [jO). 


2.2. FIR Data from AKARI and Herschel 

To obtain FIR luminosities, we first cross-identified AGNs 
against the AKARI/FIS all-sky survey bright source catalog 
(Yamamura et a l.|2010 . AKARI is the first Ja panese infrared 
astronomical satellite ( Murakami et al.|[2007) which carries 
two instruments, i.e., the Infrared Camera (I RC; |Onaka et 
al.|2007b and the Far-Infrared Surveyor (FIS;|Kawada et al. 
2007) . 

The all-sky survey has been performed with FIS in four 
bands, respectively, centered at 65, 90, 140, and 160 pm. In 
this study, we utilised 90 pm sources with the flux quality flag 
of FQUAL = 3 (i.e., high quality). The Sir-detection limit of 
the 90 pm band is 0.55 Jy. By matching the SDSS AGNs with 
the 90 pm sources within the maximum radius of 18", which 
corresponds roughly to the 3<r position error in the cross-scan 
direction of FIS, we obtained 678 AKARI/FIS counter parts 
of the type-2 AGNs, which is ~ 1.9% of all AGNs in the 
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sample. Note that we checked SDSS spectra of whole FIR- 
detected AGNs based on visual inspection in order to avoid 
unusable data, e.g., noisy data. 

To overcome the shallow flux limit of the FIS survey, we 
additionally used the PACS Evolutionary Probe (PEP) sur¬ 
vey data (Lutz et al."| |2011j >. PEP is a deep FIR photometric 
survey with the Photodetector Array Camera and Spectrom¬ 
eter (PACS; Poglitsch et al. ||2010j>, on board of the Herschel 
Space Observatory ( [Pilbratt et al.|2010) . The field selection of 
the PEP survey includes popular multi-wavelength fields such 
as GOODS, COSMOS, Lockman Hole, ECDFS, and EGS. 
Twelve objects in our AGN sample are located in the COS¬ 
MOS field ( [Scoville et al.pOOT) , and none of them were de¬ 
tected in the AKARI survey. By matching these 12 objects 
with the PEP 100 /im source catalog, which has 5cr-detection 
limit 0.0075 Jy, we obtained FIR counterparts for 11 objects. 

As described in Section 2.1 we divide our FIR sample 
into two subsamples of Seyfert 2s and LINER 2s, and re¬ 
move [O I]A6300-weak LINERs. First, we adopted criteria 
based on the BPT diagram using the [O III] A5007/H/3 and 
[O I]A6300/Ha flux ratios (i.e., Kewley et al. 2006[ > to sep¬ 
arate the sample into 355 Seyfert 2s and 231 LINER 2s, 
including six and three Herschel- detected objects, respec¬ 
tively. In this classification, we gave a criterion, S/N >3.0 
for the [O I]A6300 line. Then, adopting_a criterion, i.e.. 


JO I] A6300/H q < 1/ 6 (e.g., Filippenko & Terlevich 1992 
Alonso-Herrero et al.|2000| we removed 94 [O I]A6300-weak 
LINERs. 


In summary, we obtained 483 AKARI/FIS-detected objects 
and nine Herschel! PACS-detected objects, for which we in¬ 
vestigate the AGN-SF connection in next sections. Figure |T| 
presents the FIR-luminosity distribution of our sample as a 
function of redshift, clearly showing that PEP survey is al¬ 
most two orders of magnitude deeper and complementary to 
the shallow FIS survey sample. 


2.3. MIR Data from WISE 

We collected MIR data of 483 AKARI-detected objects, by 
matching them against the All WISE source catalog 4 , which is 
the most recent data release of the Wide-field Infrared Survey 
Explorer (WISE', [Wright et al.|2010| >, after combing all previ¬ 
ous data fro m the WISE cryogenic and N LO WISE ( |Mainzer| 
et al. 2011). We adopted the maximum radius of 6." 1, 6."4, 
675, and 1270, respectively for 3.4, 4.6, 12, and 22 fj ,m band 
images, accounting for the averaged point-spread functions 
(PSFs) in each band. For undetected sources, upper limits 
from the 5er sensitivity are given as 0.08, 0.11, 1, and 6 mJy 
respectively for 3.4, 4.6, 12, and 22 //m bands. In this process, 
we obtained the MIR fluxes for 99% of the AKARI-detected 
AGNs. 


3. RESULTS 

In this section, first, we compare four different SF indica¬ 
tors, namely, FIR luminosity, the break at 4000A (D4000), 
UV-luminosity, and the [OII] A3727 emission line luminosity, 
to investigate the reliability of the FIR luminosity as a SF in¬ 
dicator (Section [3.1| >. Second, we investigate the relation be¬ 
tween AGN and SF luminosities of the FIR-matched type-2 
AGN sample (Section [L2) . 


3.1. SF Indicators 

To test whether the FIR luminosity Lfir is a reasonable SF 
indicator, we compare Lfir to other SF indicators, i.e., D4000 
( jBrinchmann et al.|2004| ), UV luminosities, and [O II]A3727 
line luminosity. We obtained these measurements and then 
converted them to star formation rates (SFRs) as explained 
below. 

For FIR luminosity, we collected 90 and 100 /tm data re¬ 
spectively from the AKARI and Herschel samples, as de¬ 
scribed in Section [272] For given spectral energy distributions 
(SEDs) of typical SF galaxies (e.g.. Dale & Helou|2002) , the 
difference between fluxes at 90 and 100 /im is relatively small, 
e.g., log(90L’ 9 o /im )/log(100L'iooi I m) ~ 1.005. Thus, we used 
90 and 100 /tm fluxes for calculating L FIR . We also ignore 
the redshift effects since the redshift range of our sample is 
small (0.01 < z < 0.22): even for the highest redshift objects 
at z = 0.22 in our sample, the flux correction is negligible, i.e., 
log(90L 9 o Mm , obs )/ log(90 L 90Mrn , rest ) ~ 0.987. Here, we adopt 
the conversion recipe in Kennicutt (1998!: 


SFRfir (Mq year” 1 ) = 4.5 x 1(T 44 Lfir (erg s” 1 ). (1) 

Second, we obtained the SFR determined from the break at 
4000A SFR D 4 ooo from the MPA-JHU SDSS DR7 g alaxy cata¬ 
log, which is based on the technique discussed by |Brinch mann 
et al. ( |2004) . First, they constructed the relation between the 
specific SFR measured from the Ha line and D4000 using a 
sample of star-forming galaxies. Adopting this relation along 
with stellar masses estimated from the mass-to-light ratios, 
they derived SFR from D4000 for galaxies, of which emis¬ 
sion lines are not reliable as SF indicators due to the con¬ 
tamination of AGN. Note th at ap erture corrections have been 
applied ( Brinchmann et al.|2004) l, using the resolved color in¬ 
formation available for each galaxy. Since SFRd 4 ooo can be 
determined for AGN host galaxies, it has been adopted in var¬ 
ious studies (e.g.,|Netzer |2009| ). 

Third, we used the UV luminosity as a SF indicator. From 
the GALEX UV luminosities, we calculated SFRuv using the 
recipe given by Kennicutt ( |1998[ >: 

SFR UV (M q year -1 ) = 1.4 x l(T 28 Luv (erg s” 1 Hz -1 ), (2) 


2.4. UV Data from GALEX 

We also collected near-ultraviolet (NUV) and far-ultraviolet 
(FUV) flux data obtained by the Galaxy Evolution Explorer 
All-Sky-Imaging Survey (GALEX! AIS; |Martin et al. 2005). 


Using the GALEX/ AIS-SDSS matched catalogs in Bianchi et 
[ah ( 2011) , which was constructed by matching the GALEX 
GR5 data against the SDSS DR7 with a radius of 3", we ob¬ 
tained 242 and 149 counterparts out of 483 AGNs, respec¬ 
tively, in the NUV and FUV. For the remaining objects, an 
upper limit is given based on typical depths of 20.8 and 19.9 
AB magnitude in the NUV and FUV. 


where Luv is the luminosity density integrated over the spec¬ 
tral range 1500-2800A. In this section, we focus on NUV data 
because all FUV-detected objects are detected with NUV, and 
converted NUV luminosities to SFRuv- We corrected for the 
UV extinction using the Baimer decrement 5 . 

Last, we adopted the [O II] A3727 line luminosity as a SF 
indicator using the following equation (Kennicutt 11998) : 

SFRjoii] (Mq year- 1 ) = 1.4 x 10- 4I L [OU] (erg s" 1 ), (3) 

where L[oii] is the luminosity of the [O II]A3727 line. We 
corrected for the dust extinction (see footnote 5). Note that 


4 http://wise2.ipac.caltech.edu/docs/release/allwise/ 


5 http: //ned. ipac .caltech. edu/level5/S eptO 1 /Rosa/Rosa_appendix. html 
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Fig. 2.— Comparisons of the four independently estimated SFRs from FIR luminosity, D4000, UV-luminosity, and [O II]A3727 luminosity for the AKARI- 
detected (small filled and open circles) and Herschel-detected objects (stars). In panel (a), blue and red symbols represent low-z (0.05 < z) and high-z (0.05 > z) 
objects, respectively, and open symbols indicate red or old objects, i.e., D4000 > 1.8. In panels (c), (d), and (e), upper limits of the UV-based SFR are shown 
as open symbols. Gray symbols in panels (d) and (f) represent the red or old population same as open symbols in panel (a). Green-large circles mark highly 
obscured objects, i.e., Eb-v > 1. Gray-dotted line indicates the one-to-one relation, and black solid and dashed lines are fitting results, i.e., y = ax+b andy = x+c, 
respectively. Gray solid and dashed lines in panel (a) are best-fit for data excluding red or old objects. Yellow circles represent normal star-forming galaxies 
detected with AKARI. 


here we focus on S/N > 6 objects for the [O II] A3727 line 
flux to increase reliability. 

We compare the four independently estimated SFRs in Fig¬ 
ure [2] and discuss the details as follows. First, in compar¬ 
ing L FIR -based and D4000-based SFRs (Figure [2j a)), we find 
a relatively good trend at high SFR regime, while many ob¬ 
jects are below the one-to-one relation (dotted line), indicat¬ 
ing that D4000 underestimates SFR compared FIR luminos¬ 
ity. Particularly at log SFR fir < 0, the discrepancy becomes 
unacceptably large. We fit the data with a linear function, 
i.e., y = ax + b (black-solid line) or with a fixed slope, i.e., 
y = x + c (black-dashed line) as listed in Table 1. Since the 
SFR d40 oo is calibrated using starburst galaxies, the D4000 
method is subject to large uncertainties for older and red¬ 
der galaxies (see also Netzer] [2009) l. In fact, such galaxies 
(i.e., D4000 > 1.8, open symbols in Figure [2]) are mostly lo¬ 
cated below the one-to-one line. Thus, we performed a liner 
fit after excluding such galaxies (gray-solid and gray-dashed 
lines), which slightly improves the relation. Although aper¬ 
ture correction was adopted in [Brinchmann et ah (|2004| >, we 
further consider the aperture effect that the fixed 3" fiber size 
of the SDSS spectroscopy covers a smaller physical area of 
the lower-z galaxies, hence, the SFR D 4 ooo may be more under¬ 
estimated than higher-z galaxies. To test the aperture effect, 
we divided the sample into two redshift bins, i.e., z > 0.05 
(red) and z < 0.05 (blue) in Figure[2ja). However, we find no 
significant difference between them, concluding that the aper¬ 
ture effect is not the origin of the discrepancy between SFR hr 
and SFR d4 ooo- Moreover, to examine an extinction effect, we 
marked highly obscured objects, i.e., E B -v > 1 (green-large 
circles), but we found no significant extinction-related bias. 


Note that at low FIR luminosities (~ 10 43 ), the FIR is be¬ 
lieved to come from diffuse dust cirrus warmed by the back¬ 
ground starlight, not necessarily recently formed stars. In this 
case, FIR luminosity overestimates SFR although such low- 
Lpir objects are negligible in our samples (see Figure [lj. To 
demonstrate the correlation between D4000 and L F i R , we also 
plot normal star-forming galaxies detected with AKARI (yel¬ 
low circles). In conclusion, we suggest that SFRd4ooo indi¬ 
cator seems to have various issues, especially for older and 
redder galaxy populations. 


Second, we compared FIR-based and [O II]A3727-based 
SFRs in Figure [2] b). SFRs are measured and calibrated based 
on Ha and [OIIJA3727 lines, mainly for normal galaxies (e.g., 


Kennicutt 1983, 1992; Kennicutt et al. 1994} 

Madau et al. 

1998 , Hopkins et al. 2003; Kewley et al. 2004 

Moustakas et 

al. 2006). Recently, to derive attenuation-corrected line lumi- 


nosities of galaxies, some studies have combined optical and 
infrared observations (e.g., |Kennicutt et al. 2009, Dominguez 
|Sanchez et al.|[20l2| and references therein). Unfortunately, 
because of AGN contributions, it is difficult to adopt such 
corrections for our AGN sample. The best-fit relation be¬ 
tween SFRpi R and SFRjo nj (black-solid line) is steeper than 
the one-to-one correspondence, plausibly due to the combi¬ 
nation of the following effects: the AGN contribution to the 
[O II] A3727 line, leading to an overestimation of the SFR, 
particularly for high [O II]A3727 luminosity objects, and the 
stronger aperture effect of the SDSS spectroscopy for lower-z, 
i.e., lower [O II]A3727 luminosity objects. Furthermore, for 
highly reddened objects, [OII] A3727-based SFR may become 
too high if extinction is over-corrected. It would be necessary 
to consider also metallicity and ionization conditions for the 
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Fig. 3.— Relations between FIR luminosity at 90 or 100 /im, and AGN luminosity estimated from the [O III]A5007 line (the left-hand panel (a)) and the 
combination of the [O III]A5007 and [O I]A6300 lines (the right-hand panel (b)). The AKARI-detected and Herschel-detected objects are denoted with small 
circles and stars, respectively. Vertical and horizontal bars show la errors in their luminosities. Red and blue symbols indicate Seyfert 2s and LINER 2s, 
respectively. Blue-, red-, and purple-dashed lines are respective fitting results of Seyfert 2s, LINER 2s, and total objects. The reference line fro in |Nctzcr| j2009|i 
is re prese nted by black lines, assuming three different flux ratios (i.e., F60nm/Fioo Mm ). namely mean (the solid line), minimum and maximum ratios (dotted lines) 
from|Dale & Helou||2002|). The pure-AGN sequence with the la range is calculated from an intrinsic-AGN SED, shown as gray lines. Six pure-AGN candidates 
are denoted with large-ETack circles in the panel (b). Light-green circles indicate composite objects selected with the BPT diagram, which are detected with 
A KARL 


[O II] A3727 line. On the other hand, the FIR luminosity is 
contributed by dust in cirrus clouds that are warmed by dif¬ 
fuse starlight, particularly at low SFR (i.e., below 10 M ( . Jyr). 
This overestimates of SFRfir results in a steeper slope both 
for AGNs and normal star-forming galaxies. 

Third, the UV-based and FIR-based SFRs are compared in 
Figure |2}c). Many objects are located below the one-to-one 
line, indicating that UV-based SFR is largely underestimated 
presumably due to the dust extinction. Note that since UV de¬ 
tection is not available for all AGNs, we include upper limits 
of UV-based SFR (open symbols) while most obscured galax¬ 
ies are marked with green circles. The large scatter between 
UV-based and FIR-based SFRs (see Table 1), suggests that 
the UV-based SF are highly uncertain at any luminosity range, 
due to extinction and the contamination from AGB stars. 

Fourth, Figure |2jd) presents a comparison of the UV-based 
SFR with SFR D 4 ooo, illustrating no significant correlation be¬ 
tween them. On the other hand. Figure [2]e) shows a good 
correlation between the SFRuv and SFR[oiij- However, both 
SFRs suffer an extinction effect although there is a correlation 
between them. 

Fifth, Figure |2jf) compares SFRimooo and SFR[onp show¬ 
ing a correlation with a systematic shift toward high SFR[o nj. 
The larger SFR[on] than SFR D 4 () oo is mainly due to the 
AGN contribution to the [O II]A3727 line, while normal star¬ 
forming spirals show no excess of SFR[oiij. 

For quantitative analysis, we calculated the Spearman rank- 
order correlation coefficient p and their statistical significance 
p for all available data excluding upper limits. As presented 
in Table 1, we confirmed that the UV-based SFR seems to 
show week or almost no correlation with other indicators. 
FIR-based SFR seems to present slightly stronger correlations 
with optical-based SFRs than the relation among other SFRs, 
suggesting that FIR-based SFR is a reasonable SF indicator. 

By comparing four independently estimated SFRs, we con¬ 
clude that the FIR luminosity is the most reasonable SF indi¬ 
cator for host galaxies of type-2 AGNs. Thus, we will use the 


FIR luminosity as a SFR indicator, and compare it with the 
AGN luminosity in the following sections. 


3.2. A Relation between FIR and AGN Luminosities 


To examine the Lfir-Tagn relation, we need to estimate 
AGN luminosities. X-ray luminosity is the most reliable indi¬ 
cator of AGN bolometric luminosity. However, it is not avail¬ 
able for large samples including low-luminosity AGNs. The 
[Ne V]A3426 line is also a good indicator to estimate AGN 
luminosities without a contamination from H II regions (e.g., 
Gilli et al. 2010| >. Unfortunately, this emission line is typically 


weak and not covered by SDSS wavelength range for objects 
at z < 0.1. As a number of previous studies of type-2 AGNs 


luminosities (e.g., Heckman et al. 

2004 Netzer et al. 2006 

Kewley et al. 2006;]Netzer 2009 

Kauffmann & Heckman 

2009; [LaMassa et al. 2013), we calculate bolometric lumi- 


nosity from the [O III] A5007 line, adopting a bolometric cor¬ 
rection, BC = 600 (e.g., |Kauffmann & Heckman|2009 , iNetzer 
|2009| . Figure |3ja) shows the relation between the FIR lumi¬ 
nosity and AGN luminosity based on the [O III]A5007 line. 

On the other hand, Netzer ( 2009) claimed that using the 
[O III] A5007 luminosity as a proxy for the AGN bolometric 
luminosity is unreliable due to its dependence on the ioni¬ 
sation parameter, which is critical for low-ionization sources 
such as LINERs. Thus, in order to avoid this ionization effect, 
we also calculated AGN bolometric luminosity from the com¬ 
bination of [O III]A5007 and [O I]A6300 line fluxes, using the 
calibration given by Netzer (2009) : 

logTAGN = 3.53 + 0.251ogL[om]+0.751ogL[oi], (4) 


where L[oni] and T[oi] are extinction-corrected luminosities 
of [O III]A5007 and [O I]A6300 lines, respectively, in units of 
erg s -1 , using the Balmer decrement. Although the [OI] A6300 
line is much weaker than [O III] A5007, we determined reliable 
AGN bolometric luminosities for 492 objects with S/N >3.0 
for the [O I]A6300 line. Note that even in the extreme case. 
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e.g., S/N[om] = 6 and S/N[oi] = 3, the uncertainty of logarith¬ 
mic AGN luminosity, log /-agn. is 0.25, and this is acceptable 
in our discussion. The relation between the FIR luminosity 
and AGN luminosity estimated by using the combination of 
two oxygen lines in Figure [3jb). 

By comparing Figures [3ja) and |3]b), we examine which 
AGN luminosity estimates is more reliable in this study. To 
investigate the ionization parameter effect, we plotted Seyfert 
2s and LINER 2s with different symbols (i.e., red and blue 
symbols, respectively). We confirmed that the relation with 
the [O III] A5007-based AGN luminosity shows a larger scatter 
than that with AGN luminosity based on the combination of 
the [O III] A5007 and [O I]A6300 lines, particularly for LINER 
2s at low Lagn range. We performed a liner fit to Seyfert 2s, 
LINER 2s, and total objects, finding that the [O III]A5007 and 
[O I]A6300 combined method seems to correct for systematic 
trend in the distribution due to the ionization condition. By 
applying the Spearman rank-order test (see Table 1), we find 
a stronger correlation of FIR luminosity with AGN luminos¬ 
ity based on [O III]A5007 and [O I]A6300 method than based 
on [O III]A5007 only. These results imply that an ionization 
mechanism of Seyfert 2s and LINER 2s is similar. Thus, in 
the following analysis, we use both Seyfert 2s and LINER 
2s and adopt the AGN bolometric luminosity estimated by 
[O III] A5007 and [OI] A6300 lines. 


As shown in Figure [3] AKARI-detected objects show a 
clear trend between FIR and AGN luminosities, confirming 
the AGN-SF relation reported by previous studies. To di¬ 
rectly compare with the previous studies, we included the ref¬ 
erence line from |Netzer| ( |2009| ), which represents the relation 
between the D4000-based SF luminosity and AGN luminos¬ 
ity. Note that the reference line was revised by converting 
60 pm luminosity to 90 pm luminosity using three differ¬ 
ent flux ratios, i.e., \og{F a)jim /F miJXn ) = -0.55, -0.32, and 
0.21, which are based on the different SED templates (Dale 


& Helou 2002). Our result is consistent with that of |Netzer 


( 20091 within the uncertainties of the FIR luminosity conver¬ 
sion. In contrast, we do not find a strong evidence of the en¬ 


hanced SF for given AGN luminosity, as reported by Rosario 
et al. ( '2012j l for low-luminosity AGNs, particularly at high 
redshift, suggesting that low-luminosity AGNs are hosted by 
low-SF galaxies in the present day. 

In addition, we plotted Herschel- detected objects as red 
stars. Since the flux limit of these objects are two orders of 
magnitude deeper than the AKARI/FIS survey, the additional 
Herschel sample helps us to overcome the flux limit of the 
shallow AKARI/FIS survey (see Figure|T|. Herschel- detected 
objects are slightly shifted to lower Lfir/Lagn ratio compared 


to AKARI-detected objects. Note that as we described in Sec¬ 
tion [2] all type-2 AGNs in the COSMOS field are undetected 
with AKARI while most of them are detected with Herschel. 
On the other hand, ~ 30,000 objects in the SDSS field are 
not detected with AKARI, implying that a large number of 
AGNs that are not detected with AKARI, may occupy the re¬ 
gion where Herschel- detected objects are located in Figure[3] 


4. DISCUSSION 

4.1. Type-1 AGNs versus Type-2 AGNs 

We investigated the relation between FIR and AGN lumi¬ 
nosities using type-2 AGNs, for which AGN bolometric lumi¬ 
nosity is somewhat uncertain compared to type-1 AGNs. To 
overcome the uncertainty of the bolometric luminosity and to 
test whether type-1 AGNs also follow the same relation be¬ 



FlG. 4.— Relation between FIR and AGN luminosities of type-1 and type-2 
AGNs at 0.01 < z < 0.22. The gray circles, stars, and lines are same as in 
Figure[3]b). The X-ray type-2 and type-1 AGNs are represented by red and 
blue symbols, respectively. For each survey, individual symbols are given as 
labeled at the top left. In addition, SDSS type-1 AGNs with AKARI detec¬ 
tions are also shown as yellow circles. Respective fitting results are shown as 
dashed lines with same colors of each symbol. 


tween FIR and AGN luminosities, we used X-ray AGN sam¬ 
ples in this section. 

First, we collected X-ray detected type-2 AGNs at 0.01 < 
z < 0.22 from Lusso et al. ( (2011| , and matched them against 
the PEP catalog, finding six Herschel- detected sources. For 
them, we adopted AGN bolometric luminosities estimated us¬ 
ing infrared and X-ray luminosities (Lusso et a l.|2011 1 . Sec¬ 
ond, we collected X-ray detected type-1 AGNs at the same 
redshift range from the COSMOS ( |Brusa et al. 2010| l. By 
matching them against the PEP catalog, we obtained three 
Herschel- detected X-ray AGNs. Third, we collected X-ray 
detected AGNs from the 70 months Swift-BAT all-sky hard 
X-ray survey catalog ( jStem & Laor 20131. By matching them 
against the AKARI/FIS catalog, we obtained two X-ray de¬ 
tected type-1 AGNs and four X-ray detected type-2 AGNs. 
For these AGNs, we estimated AGN bolometric luminosity 
from the X-ray luminosity with a bolometric correction in 
|Rigby et~ak ( |2009| l. Fourth, we collected 23 Sey fert Is and 
18 Seyfert 2s from the 12 micron galaxy sample ( |Spinoglio| 
|& Malkan|1989| [Rush et al.|1993[), which are d etected in X- 
ray with XMM-Newton ( |Bri g htman & Nandra]|201 1| >. Using 
these four samples, we plotted 28 type-1 AGNs (blue) and 28 
type-2 AGNs (red) detected in X-rays on Figure [4] These X- 
ray AGNs seem to generally follow the similar trend between 
FIR and AGN luminosities, albeit the small sample size. To 
quantitatively assess these trends, we plotted the best-fit lin¬ 
ear relations, respectively for X-ray type-2 AGNs and X-ray 
type-1 AGNs, and calculated the Spearman rank-order corre¬ 
lation coefficients (see Table 1). We find that X-ray type-2 
AGNs show the similar relation to our SDSS type-2 AGNs, 
while X-ray type-1 AGNs seems to have a different slope. 
The shallow slope of X-ray type-1 AGNs seems due to the 
narrow Lagn range. A larger sample is required to examine 
the origin of the difference in slope. Note that systematic er¬ 
rors from different data sets, especially in estimations of AGN 
luminosities, would affect these comparisons. 

In addition to X-ray AGNs, we selected optically-selected 
type-1 AGNs using the SDSS low-luminosity AGN sample in 
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Fig. 5.— Redshift distributions of our type-2 AGNs on the Lfir-^agn 
plane. The sample is divided in three redshift bins, 0.01 < z < 0.04, 0.04 < 
z < 0.10, and 0.10 < z < 0.22, respectively denoted with black, blue, and red 
symbols. The AKARI-detected and Herschel -detected objects are shown as 
open circles and filled stars, respectively. The luminosity limits at z = 0.01, 
0.04, and 0.10 based on the AKARI 5cr-detection limits are denoted with 
horizontal dashed lines. Gray lines are the same as those in Figure [3] 


|Stern & Laor| ( |2013] ), for which the broad Ha line is detected, 
enabling us to estimate AGN luminosity b ased on the H a line 
luminosity based on a recipe in Greene & Ho| (j2005 i. By 


matching them against the AKARI/FIS catalog, we obtained 
45 type-1 AGNs. For these objects, an AGN bolometric lu¬ 
minosity is estimated from the broad Ha line luminosity. In 
Figure [4] we plotted optical type-1 AGNs along with type-2 
AGNs with yellow circles, showing that optical type-1 AGNs 
follow a consistent relation between FIR and AGN luminosi¬ 
ties. We conclude the method of AGN bolometric luminosity 
estimation, i.e., narrow emission line luminosity. X-ray lumi¬ 
nosity, and broad Ha luminosity, does not significantly affect 
the relation between FIR and AGN luminosities, and that op¬ 
tical type-1 and type-2 AGNs show a similar relation between 
FIR and AGN luminosities. 

4.2. Comparisons with the Previous Studies 

In this section we compare our result with that of the pre¬ 
vious studies. As shown in Figure[3j we found a similar rela¬ 


tion between SF and AGN luminosities as reported by Netzer 
(20091. The difference between our study and that of Netzer 
(2009) is that while we used FIR luminosity as a proxy for 


SF, Netzer ( |2009 1 used D4000 in estimating SF luminosity. 
As discussed in Section [XT] D4000-based SFR is not well de¬ 
termined at lower SFR since the calibration was based on star- 
burst galaxies. Even with the FIR luminosity, which is a rel¬ 
atively better SF indicator, we find a similar relation between 
SF luminosity and AGN luminosity. This is probably due to 
the fact that the relation is not tight and the ratio between 
SF luminosity and AGN luminosity has a broad distribution, 
hence, the systematic difference between D4000-based SF lu¬ 
minosity and FIR luminosity is not clearly detected. 

It is interesting to note that Herschel -detected objects fol¬ 
low the similar relation as AKARI-detected sources, i.e., low- 
L A gn AGNs at higher redshift, show the similar trend between 
FIR and AGN luminosities, suggesting that the relation is not 
due to the selection effect. In Figure [5] the luminosity lim¬ 
its with increasing redshift is demonstrated in the Lfir-Lagn 
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Fig. 6.— Mean FIR luminosities of the AKARI-detected objects for each 
L A gn bin- The mean luminosities in each redshift range, 0.0l < z < 0.04, 
0.04 < z < 0.10. and 0.10 < z < 0.22, are denoted with black, blue, and red 
circles, respectively. For each mean value, vertical bars represent 3a errors 
while horizontal bars show the AGN-luminosity ranges. If the sample size is 
less than 2 in a bin, a filled circle without error bars is given. The Herschel- 
detected objects are shown as filled stars. Horizontal and gray lines are the 
same as in Figure[5] 


plane. Here, we divided our sample into three redshift bins, 
i.e., 0.01 < z < 0.04, 0.04 < z < 0.10, and 0.10 < z < 0.22, 
and the AKARI/FIS 5<r-detection limits at z = 0.01, 0.04, and 
0.10 are denoted with dashed horizontal lines. As shown in 
Figure [5] the AKARI-detected sources are strongly affected 
by the Malmquist bias, indicating that the AKARI/FIS sample 
alone does not allow us to investigate the relation between FIR 
and AGN luminosities without suffering the selection effect 
due to the flux limit. In contrast, Herschel- detected sources 
enable us to examine the relation over a wide luminosity range 
at given redshift (e.g., 0.1 < z < 0.22 or 0.04 < z < 0.1). In 
particular, for AGNs at 0.01 < z < 0.22, the relation between 
FIR and AGN luminosities is detected over a wide range of 
AGN accretion luminosity, 42 < logL A GN (erg s -1 ) < 46. We 
conclude that the relation between FIR and AGN luminosities 
is not due to the flux limit of the AKARI/FIS surveys. 

To investigate the effect of the flux limit for given redshift 
bins, we calculated the mean FIR luminosities of the AKARI- 
detected sources in each redshift bin, i.e., 0.01 < z < 0.04, 
0.04 < z < 0.10, and 0.10 < z < 0.22, after dividing the AGNs 
in each redshift bin into subgroups based on AGN luminosity. 
In Figure [6] we present the mean FIR luminosities for each 
bin. The mean FIR and AGN luminosities show rather a flat¬ 
tened pattern in each redshift bin, compared to the relation 
between FIR and AGN luminosities of individual objects. In 
particular, at low AGN luminosity, the mean SF luminosity 
appears to be enhanced for fixed AGN luminosity, as simi¬ 
lar reported by recent studies (e.g., Sha o et al.|2010| [Rosario 
et al.||2012l. However, this flattened pattern is not detected 


when we used individual luminosity measurements instead of 
mean luminosities (see also Netzer 2009) ). The reason why 
we do not find AGNs hosted by galaxies with enhanced SF 
(i.e., above the one-to-one relation in Figure [3]) may result 
from our sample selection since we excluded the composite 
objects in the BPT selection. It is possible that we missed SF- 
enhanced AGNs, which could be classified as composite ob¬ 
jects, i.e., star-burst AGNs. To check whether composite ob- 
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Fig. 7.— Composite SEDs of pure-AGN candidates (red) and star-forming AGNs (black), normalized with SDSS R-band luminosity. Three representative 
SEDs are also shown as gray lines: starburst galaxy IRAS 19254—7245 (solid line), NGC 6240 (dashed line), and an average Seyfert 2 galaxy (dash-dot line). 


jects are distributed above our relation or not, we plot galax¬ 
ies classified as composite objects based on the BPT method 
(shown as light-green circles in Figure [3). Note that AGN 
bolometric luminosities might be overestimated for compos¬ 
ite objects because of the contamination from star formation. 
As shown in Figure [4] we investigated the Lagn-^fir relation 
of X-ray selected objects and confirmed that several X-ray de¬ 
tected AGNs are distributed on the enhanced-SF area (top-left 
region in this figure), although its fraction is quite low, in¬ 
dicating that our results are consistent to Netzer ( |2009j ). A 
direct test with the SDSS composite objects is difficult to per¬ 
form since the AGN luminosity estimated from the narrow 
emission lines would be much more uncertain due to the con¬ 
tribution from SF. Note that Rosario et al. ( |2012| adopted 60 
pm luminosity as a SF indicator which may be more contam¬ 
inated by an AGN component than 90 and 100 pm luminosi¬ 
ties (e.g., |Spinoglio et al.|2002| , and their mean FIR luminosi¬ 
ties may be overestimated. Moreover, if Rosario et al.| ( [20T2| > 
lost low X-ray luminosity objects which would show low FIR 
luminosities in their sample selection, their averaged FIR lu¬ 
minosities would be overestimated although they calculated 
mean FIR luminosities by considering FIR detected and un¬ 
detected objects. 


4.3. AGN Contribution to FIR Luminosities 

If there are luminous AGNs hosted by low-SF galaxies, we 
may find them at the bottom right of the L fir -Lagn plane (see 
Figure [3}. Usually, the AGN contribution to FIR is believed 
to be negligible, since SF galaxies dominates at FIR. When 
very low FIR luminosities are probed, however, it is neces¬ 
sary to quantify the AGN contribution. To investigate the 
FIR luminosities of a pure-AGN without SF, we adopted a 
SED template from jMullaney et akj ( |201lj ). Using the strong 
correlation between the MIR (12.3 pm) and the X-ray (2-10 
keV) luminosities from |Gandhi et al.| ( |2009| ), we obtained the 
MIR luminosity as a function of AGN luminosity (see also 
|Ichikawa et al.|20f2 ), then calculate the FIR luminosity based 
on the pure-AGN SED template. In this process, we estimated 
the AGN bolometric luminosity from the X-ray luminosity 


(Rosario et al.|2012] >. 

In Figure [3] the estimated pure-AGN sequence is denoted 
with gray lines. Based on a comparison our sample with the 
pure-AGN sequence, we found that the AGN contribution in 
our sample seems to be negligible, although there are six ob¬ 
jects reaching this pure-AGN sequence. Note that it is im¬ 
portant to examine these objects located on the pure-AGN se¬ 
quence since they are likely to be low-SF AGNs compared to 
the star-forming galaxies in our sample, if the expectation of 
contribution to FIR luminosities from the intrinsic-AGN SED 
is correct (see also |Rosario et al.|2012j[Mor & Netzer|2012) >. 


4.4. Spectral Energy Distributions 

For understanding the L F ir-Tagn relation in detail, in this 
section we investigate SEDs of our type-2 AGNs. As we 
mentioned in Section |4.3| six objects are located close to 
the pure-AGN sequence. Here, we focus on these objects as 
pure-AGN candidates, i.e., luminous AGNs hosted by no- or 
low-SF galaxies (see large circles in Figure [3]>. Using multi¬ 
wavelength data, i.e., SDSS-optical, AKARI-FIR, WISE- 
MIR, and GALEX-UV data, we separately constructed the 
composite SEDs of pure-AGN candidates and AGNs hosted 
by star-forming galaxies shown in Figure[7] First, we normal¬ 
ized all-band luminosities with the SDSS A 1 -band luminosity. 
Note that this normalization is effectively a stellar mass nor¬ 
malization, since optical-NIR luminosities roughly represent 
the stellar mass. In this construction, we used 5cr upper lim¬ 
its for undetected sources. Along with the SEDs, we plotted 
three SED templates: a starburst galaxy IRAS 19254-7245, 
NGC 6240 (starburst + Seyfert 2), and a composite SED of 
Seyfert 2 galaxies (Polletta et al. 2007j l as gray lines in Fig¬ 
ure [7] These template SEDs are also normalized at the SDSS 
/x’-band wavelength. 

As shown in Figure[7] we confirmed that the FIR luminosi¬ 
ties of pure-AGN candidates are lower than AGNs hosted by 
star-forming galaxies, as already suggested in Figure [3] This 
indicates that pure-AGN candidates have significantly lower- 
SF host galaxies relative to star-forming AGNs. As SF indi¬ 
cators, UV luminosities also show a similar trend to the FIR 































FIR and Accretion Luminosities of the Present-Day AGNs 


9 



1 10 
Rest Wavelength (/xm) 


45.0 < log i AGN ,[oui] 4£ [oi] < 46 0 


Pure-AGN candidates (3) 
Star-forming AGNs (115) 



Rest Wavelength (/xm) 



0.1 1 10 100 
Rest Wavelength (/xm) 


Pure-AGN candidates (2) 
star-forming AGNs (197) 




Fig. 8.— Composite SEDs of pure-AGN candidates and star-forming AGNs, for three AGN-luminosity bins (left-hand panels), and for three FIR-luminosity 
bins (right-hand panels). The lines are the same as in Figure]?] 


luminosities, i.e., UV luminosities of pure-AGN candidates 
are lower than star-forming AGNs, although UV luminosities 
suffer from uncertainties as described in Section 13.11 More¬ 
over, in the optical range, there is a difference of the spectral 
slope between pure-AGN candidates and star-forming AGNs, 
suggesting the difference of D4000 between two groups. On 
the other hand, MIR luminosities, which are an indicator of 
the AGN luminosity, show no significant difference between 
pure-AGN candidates and star-forming AGNs. These results 
indicate that for given stellar mass, the AGN luminosity is 
comparable between pure-AGN candidates and star-forming 
AGNs while pure-AGN candidates have on average lower SF 
than star-forming AGNs. Compared to the templates, we 
found the SED of pure-AGN candidates is located between 
NGC 6240, which is a composite of starburst and Seyfert 2, 
and Seyfert 2 templates while the SED of star-forming AGNs 
is similar to the template of NGC 6240. This indicates that 
pure-AGN candidates are hosted by significantly lower-SF 
galaxies than star-forming AGNs. Note that even pure-AGN 
candidates have higher FIR luminosity than Seyfert 2 tem¬ 
plate, probably due to the AKARI flux-limit. In other words, 
AKARI sample is biased toward active star-forming hosts. 

To understand the dependency of the SEDs in relation 
with Lfir and L AG n, we divided our sample into subsam¬ 
ples by using FIR luminosities or accretion luminosities: we 
adopted AGN-luminosity bins of 43.0 < logLAGN < 44.0, 
44.0 < logL AGN < 45.0, and 45.0 < logL AGN < 46.0, and for 


FIR-luminosity bins we used 42.5 < logL FIR < 43.5, 43.5 < 
log Lfir < 44.5, and 44.5 < logLpiR < 45.5. For each bin, we 
constructed an average SED as plotted in Figure [8] 

In the left hand panels, the SEDs for three L AGN - bins are 
shown from bottom to top with increasing accretion luminosi¬ 
ties. In the case of star-forming AGNs shown as black lines, 
FIR and UV luminosities are increasing with increasing AGN 
luminosities, as expected from the positive trend between FIR 
and AGN luminosities. Also, we found that more luminous 
objects show slightly flatter slopes in the optical range, again 
suggesting the correlation between SF and AGN luminosities 
at fixed stellar mass. Moreover, we confirmed that the WISE 
continuum seems to increase with AGN luminosities, indi¬ 
cating that MIR bands would also be good L AGN indicators 
(e.g., Gandhi et al.|2009) . For each AGN luminosity bin (top 
and middle panels) the SEDs of pure-AGN candidates show 
lower-SF activities than star-forming AGNs, which is consis¬ 
tent with the result in Figure [7] 

When we divided the sample into three Lfir bins as plotted 
in the right hand panels, we found that FIR and UV luminosi¬ 
ties (i.e., SF indicators), and MIR luminosities (i.e., a tracer 
of AGN accretion) are increasing together as expected from 
the positive Lfir-L AG n trend. The spectral slopes in the op¬ 
tical range are increasing with increasing FIR luminosities, 
reflecting the decrease of the D4000 (increase of SF). Note 
that the spectral slopes in the MIR range become steeper with 
increasing FIR luminosities, implying that there is a SF con- 
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Fig. 9.— Simulations of the Lfir-Lagn distribution, demonstrating the effect of the flux limits in the AKARI/FIS all sky survey (left-hand panel), and the 
limited volume of the Herschel/PACS survey (right-hand panel). The logarithmic number density is calculated for each area bin with a size AlogLpiR = 0.2 and 
AlogLAGN =0.1, and represented with different colors. For comparison, the AKARI-detected and Herschel -detected sources are also plotted with gray circles 
and red stars, respectively, while the pure-AGN sequence is denoted with white lines. 


tribution to MIR luminosities. In the case of pure-AGN candi¬ 
dates (middle and bottom panels), pure-AGN candidates show 
higher MIR luminosities than star-forming AGNs, as expected 
from the trend in Figure [3] that pure-AGN candidates have 
higher AGN luminosities than star-forming AGNs at fixed 
FIR luminosities. 

Based on the SED analysis with FIR and AGN luminos¬ 
ity bins, we confirmed that pure-AGN candidates are hosted 
by low-SF galaxies compared to star-forming AGNs. Thus, 
these pure-AGN candidates could be a crucial sample for un¬ 
derstanding the AGN-SF connection. Since the fraction of 
such objects appears to be small, ~ 1%, this population may 
not be dominant in black hole growth history. However, it is 
possible that we are missing low-SF AGNs on the Lfir-Lagn 
plane due to the observational limitations (see Section 4.5 1 . 


4.5. The Effects of the Flux Limit and Volume Limit 

To better understand the observed relation between FIR and 
AGN luminosities in Figure [3] in this section we investigate 
the effects of the flux and volume limits by simulating the 
number density in the Lfir-Lagn plane. Observationally it is 
difficult to find objects with high Lagn and low Lfir due to 
the following two effects: 

- at lower redshift, e.g., z < 0.04, it is difficult to detect 
high-LAGN AGNs due to the limited survey volume. 


- at higher redshift, e.g., 0.10 < z < 0.22, objects with 
high Lagn are easier to detect, however, the flux lim¬ 
its of the FIR survey prevent the detection of 1ow-L F ir 
galaxies hosting high-LAGN AGNs. 


To quantify these effects, we simulate the distribution of 
AGNs as a function of AGN luminosity using the [O III] A5007 
luminosity function derived from the COSMOS and SDSS 
type-2 AGN samples (Bongiorno et al .j2010( l. First, we calcu¬ 
lated survey volumes for each redshift bin using the survey ar¬ 
eas of the SDSS DR7 (i.e., 9380 deg 2 ) and the PEP-COSMOS 
field (i.e., 2.0069 deg 2 ), respectively for the AKARI and Her¬ 
schel samples. Then, using the flux limits, i.e., 0.55 Jy for 
the AKARI/FIS survey and 7.5 mJy for the Herschel/ PACS 
survey, we estimated object numbers for each area box with 


a fixed size (AlogLpiR = 0.2 and AlogLAGN = 0.1), with in¬ 
creasing redshift. Finally, we integrated the number of objects 
over the redshift range 0.01 < z < 0.22. 

Figure |9]presents the simulation results for the AKARI sur¬ 
vey and the Herschel survey along with the observations. For 
simplicity, we adopted a [O III]A5007 line as a proxy for the 
AGN bolometric luminosity with a bolometric correction, BC 
= 600 (e.g., Kauffmann & Heckman|2009] Netzer 2009 1 , since 
we used the [O III ] A5007 luminosity function. The simulated 
distribution well reproduces the observed Lfir-Lagn relation. 
Note that since we did not include the FIR luminosity function 
in our simulations, the top-left area shows very high object 
numbers. However, if we apply a FIR luminosity function, 
the object number in this area would decrease to zero. 

Our simulation clearly indicates that the flux limit of the 
AKARI/FIS survey is insufficient to explore AGNs hosted 
by low-SF galaxies, close to the pure-AGN sequence while 
the limited volume of the PEP-COSMOS survey prevents us 
from detecting high Lagn and low Lfir sources. We conclude 
that the combination of the AKARI-detected and Herschel- 
detected AGNs used for our investigation suffers the observa¬ 
tional limitations due to the flux limit and the survey volume. 
Thus, we were not able to investigate the number density of 
the pure AGNs or AGNs in post-starburst galaxies with the 
AKARI/FIS and Herschel/ PACS data. For better understand¬ 
ing the Lfir-Lagn relation, it is required to have a wide and 
deep FIR surveys with the next-generation infrared astronomy 
missions, e.g., the Space Infrared Telescope and Cosmology 
and Astrophysics (SPICA). 


5. SUMMARY AND CONCLUSIONS 

To understand the AGN-SF connection, we investigated the 
relation between AGN and SF luminosities for a sample of 
SDSS type-2 AGNs at z < 0.22, based on the AKARI/FIS 
all-sky survey and the PEP COSMOS survey. We estimated 
AGN luminosities from [O III]A5007 and [O I]A6300 emis¬ 
sion lines, and utilized the proposed linear proportionality of 
SFR with FIR luminosities in Kennicutt’s equation. The main 
results are summarized as follows. 


1. By comparing four independent SF indicators, i.e., FIR- 
based, UV-based, D4000-based, and [O II]A3727-based 
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SFRs, we find that the FIR luminosity is the most 
acceptible and less subjective to AGN cont amination 
compared to other SF indicators (Section [3T| . 

2. There is an apparent positive trend between FIR and 
AGN luminosities for local type-2 AGNs. In contrast 
to other studies ( [Rosario et al.|2012) >, we find that low- 
Lagn AGNs also follow the similar relation between 
FIR and AGN luminosities (Section [L2] |. 

3. Using X-ray AGNs and optical type-1 AGNs, we find 
a similar relation between Lfir and Lagn, suggesting 
that the observed relation is not significantly affected 
by the method and uncertainty of the AGN bolometric 
luminosity estimation (Section |4~Tj ). 

4. The flux limit of AKARI FIR survey significantly af¬ 
fects the distribution in the L fir -Lagn plane while the 
deep-FIR data such as Herschel survey can overcome 
the limitation. It is possible that the AKARI FIR- 
detection limit is responsible for the observed trend in 
Lfir and Lagn plane (Section |4~2} . 

5. FIR luminosities of most type-2 AGNs are dominated 
by non-AGN continuum, while the AGN contributeion 
to FIR emission is negligible (Section |4.3[ i. 

6. Based on the simulation of the AGN number distribu¬ 
tion, we showed that the observed L F ir-Lagn relation 
can be explained by the flux limit of the AKARI/FIS 
survey and the limited volume of the Herschel /PACS 
survey, demonstrating the limitations of the current 
survey data for detecting and investigating luminous 
AGNs hosted by low-SF galaxies (Section [43] ). 

Although it is possible that the observational limitations may 
cause an artificial correlation between FIR and AGN lumi¬ 
nosities, the observed positive relation may suggest an intrin¬ 
sic connection between SF and AGN activities in the present- 
day, implying that the growth of stellar mass and black hole 
mass are linked at least in the AGN phase in galaxy evolu¬ 
tion. Quantifying the number density of luminous AGNs in 
low-SF or non-SF galaxies requires wide and deep future FIR 
surveys, e.g., SPICA. 
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